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SECTION  1 
INTRODUCTION 


Using  tapered  charges  is  currently  the  best  way  to  simulate 
the  pressure  pulse  generated  by  an  underwater  nuclear  detonation 
(Reference  1).  At  depths  less  than  approximately  30  m  (100  ft) 
the  simulation  is  quite  good,  if  the  pressure  pulse  is  suffi¬ 
ciently  planar  at  the  model  location,  and  if  the  explosion  bubble 
vents  through  the  surface. 

Recent  survivability  requirements  (Reference  2)  demand 
exposure  of  models  to  simulated  nuclear  pressure  pulses  at  depths 
of  300  to  600  m  (1,000  to  2,000  ft).  As  part  of  the  current 
Underwater  Explosion  (UNDEX)  Test  Series,  l/9th-scale  models  will 
be  tested  in  deep  water  to  both  conventional  high  explosive  and 
simulated-nuclear  pressure  pulses  (Reference  3). 

The  effects  of  reloading  due  to  the  repeated  explosion  bub¬ 
ble  pulses  at  these  depths  are  not  well  known,  and  are  generally 
much  more  severe  for  conventional  charges  than  for  nuclear 
detonations.  Since  most  models  to  be  tested  are  impulse 
sensitive,  understanding  the  pulsation  of  the  explosion  bubble  is 
becoming  increasingly  important,  even  for  conventional  depth 
charqes.  For  deep  underwater  nuclear  shock  simulations  using 
tapered  charges,  very  little  is  known  about  the  characteristics 
of  the  accompanying  explosion  bubble,  and  bubble  reloading  must 
be  substantially  reduced  or  eliminated  (Reference  4). 
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Physics  International  Company  (PI)  participated  in  a  re¬ 
search  effort  to  examine  explosion  bubble  growth  ana  collapse 
from  conventional  explosives  in  deep  water.  Calculations  were 
performed  for  both  spherical  and  tapered  high  explosive  charges, 
and  methods  were  examined  by  which  bubble  pulsation  might  be 
eliminated,  or  at  least  substantially  reduced.  This  report 
represents  an  initial  effort  to  calculate  and  assess  the  effects 
of  nonspherical  bubble  pulsation  in  deep  water,  a  subject  that 
may  be  of  increased  interest  during  the  1980s. 

In  Section  2  we  review  tapered  charge  designs,  and  present 
the  "baseline"  tapered  charge  geometry  used  in  the  two-dimen¬ 
sional  calculations;  in  Section  3  we  discuss  the  equation-of- 
state  models.  In  Section  4  we  review  the  literature  concerning 
explosion  bubbles  from  spherical  explosive  charges,  and  verify 
that  bubble  growth  and  collapse  over  the  first  bubble  period  can 
be  accurately  computed  for  this  simple  case.  In  Section  4  we 
also  review  what  is  known  about  nonmigrating,  multiple  bubble 
pulsation.  In  Section  5  we  present  the  calculation  of  the 
tapered  charge  explosion  bubble,  and  conclude  that  bubble 
pulsation  will  occur.  Three  methods  of  bubble  pulse  mitigation 
are  then  examined.  In  Section  6  we  present  the  results  of 
investigating  the  originally  proposed  method,  that  of  supplying  a 
propellant  gas  source  to  the  explosion  bubble  to  prevent,  or  at 
least  cushion  bubble  collapse.  In  Section  7  we  discuss  two 
additional,  more  promising  methods  of  bubble  pulse  mitigation. 
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SECTION  2 


BASELINE  TAPERED  CHARGE  DESIGN 


Tapered  charges  have  been  used  for  more  than  20  years  to 
simulate  nuclear  shock  waves  in  water.  A  typical  nuclear  shock 
wave  in  water  is  shown  in  Figure  2.1  (Reference  5).  Tapered 
charges  are  long,  slender  explosive  charges;  Figure  2.2  (Refer¬ 
ence  5)  shows  an  example  of  such  a  charge.  The  charge  is  deto¬ 
nated  at  one  end,  and  the  time  for  complete  detonation  of  the 
charge  is  controlled  by  the  total  charge  length.  Figure  2.3 
shows  water  pressure  and  velocity  for  one  such  charge  detonated 
in  shallow  water  (Reference  5).  The  gages  were  located  along  the 
axis  of  the  charge,  off  the  charge  end  where  the  charge  was  de¬ 
tonated.  A  plateau  pressure  of  2  MPa  (300  psi)  was  obtained; 
this  pressure  can  be  adjusted  by  varying  the  total  explosive 
weight  and/or  the  distance  of  the  gage  location  from  the  charge 
(Reference  6). 

Tapered  charges  have  been  built  and  expended  with  lengths 
from  0.6  to  10.7  m  (2  to  35  ft)  and  total  explosive  weights  of  4 
to  4,500  kg  (10  to  10,000  lb).  The  most  common  types  have 
weights  of  227  to  454  kg  (500  to  1,000  lb),  diameters  of  0.3  to 
0.46  m  (12  to  18  in.),  and  lengths  of  approximately  6. 1  m  (20  ft) 
(Reference  7).  An  up-to-date  listing  of  information  on  tapered 
charges  fired  in  the  past  is  maintained  by  Weidlinger  Associates, 
Chesapeake,  VA  (Reference  4). 
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Figure  2.1  Measured  nuclear  underwater  shock  wave  (courtesy 
of  Mr,  J.  D.  Gordon,  UERD) . 
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Figure  2.3  Water  pressure  and  velocity  from  tapered  charge 
detonated  in  shallow  water  (courtesy  of  Mr, 

John  D.  Gordon,  UERD) . 
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A  typical  design  was  obtained  from  personnel  of  the  Under¬ 
water  Explosions  Research  Division  ( UERD)  of  the  David  W.  Taylor 
Naval  Ship  Research  and  Development  Center,  Portsmouth,  VA  (Ref¬ 
erence  7).  This  design  is  shown  in  Figure  2.4;  this  charge  was 
modeled  in  the  two-dimensional  (2D)  calculations.  A  nominal 
water  depth  of  305  m  (1,000  ft)  was  chosen  for  the  depth  of  the 
detonation. 
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4  Representative  design  for  a  454-kg  (1000-lb) 

tapered  charge  (courtesy  of  Mr.  J.  Krezel,  UERD) 


SECTION  3 


EQUATIONS  OF  STATE 


To  perform  the  calculations  of  bubble  growth  and  shock  wave 
propogation  resulting  from  the  detonation  of  a  tapered  charge 
under  water,  equations  of  state  (EOSs)  for  High  Blast  Explosive-1 
( HBX-1 )  and  water  were  required.  An  EOS  for  TNT  was  also  re¬ 
quired  for  the  ID  check  problems  performed  to  compare  the  bubble 
growth  and  shock  wave  parameters  with  the  empirically  predicted 
values.  Because  the  existing  HBX-1  EOS  was  found  to  be  inade¬ 
quate  for  underwater  explosion  bubble  problems,  the  TNT  EOS  was 
used  for  the  2D  calculational  effort  as  well.  This  section  pre¬ 
sents  the  equations  of  state  used  for  these  materials. 

3.1  EXPLOSIVES 


The  EOS  for  TNT  is  well  described  by  the  Jones-Wilkins-Lee 
( JWL)  model  (Reference  8).  This  EOS  describes  the  adiabatic 
expansion  of  the  explosion  products  for  a  wide  range  of  explo¬ 
sives,  including  TNT.  The  JWL  EOS  form  is  widely  accepted,  and 
will  not  be  described  again  here,  as  it  is  very  well  described  in 
Reference  8.  The  TNT  coefficients  used  are  given  in  Table  3.1. 

HBX-1  is  an  aluminized  explosive,  with  constituents  as  shown 
in  Table  3.2  (Reference  9).  An  EOS  was  developed  for  this  explo¬ 
sive  by  Roslund  and  Coleburn  of  the  U.S.  Naval  Ordinance  Labora¬ 
tory  (Reference  10).  The  EOS  form  is  a  modified  JWL  form  where 


spam 


Table  3.1  JWL  EOS  coefficients  for  TNT. 

A  =  3.712 

B  =  0.0323 

C  =  0.0104527 

Rx  =  4.15 

R2  =  0.95 

b)  =  0.30 

E0  =  4.46  X  10 3  J/g  (0.07  Mbar-cm3/cm3 ) 

pQ  =  1.56  Mq/m3  (1.56  g/cm3) 

D  =  6.93  m/ms  (0.693  cn.-'ys) 

Eq:  Total  energy  released  per  unit  explosive  mass 
pQs  Initial  explosive  density 
D:  Explosive  detonation  velocity 
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Table  3.2  Constituents  of  HBX-1 


Constituent  Percent 

RDX  40 

TNT  38 

Aluminum  powder  17 

D-2  wax  5 

Calcium  chloride  0.5 


13 


the  pressure,  P,  in  the  explosive  is  given  as  a  function  specific 
volume,  V: 

(B  exp  (-kV)  +  WG.V”(1+W)  V  <  V 

P  =  _i  C  (3.1) 

(g2  (V  +  V*)  V  >  vc 

The  constants  of  Equation  3.1  and  Chapman- Jouguet  (C-J)  proper¬ 
ties  for  HBX-1  are  given  in  Table  3.3.  These  explosive  proper¬ 
ties  were  determined  largely  through  tests  in  which  the  explosive 
products  expanded  into  air.  Figure  3.1  compares  the  TNT  and 
HBX-1  adiabats.  Because  of  the  late-time  contribution  of  the 
aluminum  powder  to  the  explosive  energy  release,  the  pressures  at 
large  relative  volumes  (greater  than  100)  are  about  seven  times 
larger  for  HBX-1  than  for  TNT.  It  is  the  pressures  at  these 
large  specific  volumes  that  control  the  bubble  growth  at  time? 
when  the  bubble  radius  is  close  to  the  maximum  radius. 

3.2  WATER 

The  water  EOS  chosen  was  that  due  to  Walker  and  Sternberg 
(Reference  11).  This  model  gives  pressure  as  a  function  of  den¬ 
sity  and  internal  energy  on  the  Hugoniot,  and  for  adiabatic  re¬ 
lease  off  the  Hugoniot  for  initial  shock  pressures  up  to  25  GPa 
(250  kbar).  Table  3.4  summarizes  the  EOS  in  code  units.  Isen- 
tropes  for  this  EOS  are  compared  with  those  obtained  from 
Butkovich  (Reference  12)  in  Figure  3.2.  The  Walker  and  Sternberg 
EOS  makes  no  attempt  to  model  the  production  of  steam  below  the 
saturation  line,  so  the  two  models  disagree  below  that  point. 
Above  the  saturation  line,  however,  the  two  models  agree  very 
well.  The  question  of  steam  production  is  discussed  in 
Section  4. 
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Table  3.3  EOS  constants  for  HBX-1 


p  =  1.712  g/cm3 
pQ  =  1.624  g/cm3 
D  =  7.307  mm/ps 


YCJ 

PCJ 

PCJ 

PCJ 
dD/dp  D 

a 

W 

VC 

PC 

k 

B 

G1 

G2 

V* 


2.934 

22.04  GPa  (220.4  kbar ) 

1.858  mm/ys 

2.178  g/cm3 

3800  m/s/g/cm3 

0.1325 

0.3432 

1 . 15cm3/g 

0.3645  GPa  (3.645  kbar) 
6.450  g/cm3 
418.7  GPa  (4.187  Mbar) 
0.397  GPa  (3.97  kbar) 
9.142  x  10^  (cm/s)^ 
1.358  cm3/g 


RELATIVE  VOLUME,  V  (pQ/p) 


Figure  3.1  JWL  TNT  adiafcat  compared  with  the 
modified  JWL  adiabat  for  HBX-1. 


Table  3.4  Summary  of  water  EOS  model  of  Walker  and  Sternberg 


p  (mbar)  ■  p  +  f2  p3  +  fj  p5  +  p7 
p  »  current  density  (g/cm3)(  p0  »  0.9982  g/cm3 


0  *  E  i  0.006 


0.006  <  E  <  0.017 
<  E'  -  E  -  0.006) 


0.017  <  E 

(E*  -  E  -  0.017) 


0  <  E  <  0.0032 


0.0032  <  E  <  0.0245 
( E 1  -  E  -  C.0032) 


0.0245  <  E 

(E“  »  E  -  0.0245) 


»  0.005722427  -  1.240522  E  +  50.42535  E2 

-  1.400579  X  103E3  +  4.137950  x  106E4 

-  2.726437  x  10SE5  -  1.295684  X  10UE6  +  1.437988  x  1013E7 

fj  -  0.001015091  -  0.3270122  E’  +  6.734616  E'2 

♦  1.552785  X  104E'3  -  2.926440  x  106E'4  +  2.139341  x  108E'5 

-  ‘>.61538  x  loV6 

■  ‘>.607572  x  10“4  +  0.1122840  E* 

+  >.27579  E"2  +  82.21745  E*3  -  147.1514  E’4 

-  1.044093  x  103E*5  -  3.130131  x  104E*6 

t2  *  '  0.02748180  +  1.691130  E  +  17.12981  E2 

♦  1 . <?3364  x  104E3  -  1.549072  x  107E4 

♦  3.415591  X  109E5  -  2.357818  x  10UE6 

f2  -  -  0.02215430  +  1.510990  E’  -  10.56299E'2 

-  5.411856  x  103E'3  +  6.176871  x  105E’4 

-  1.810118  x  107E*5  -  6.205700  x  108E'6  +  4.406075  x  1010E'7 

-  6.587460  X  10UE’8 

£2  •  0.002499950  +  0.9374720  E“ 

-  4.624610  E*2  -  44.52203  E*3  +  375.1364  E*4 

f3  «  0.0268  -  0.4148  E 
f4  -  -  0.005  +  0.0741  E 


E  -  internal  energy  in  mbar-cm3/g 


23 


SECTION  4 


UNDERWATER  EXPLOSION  BUBBLE  GROWTH  AND 
COLLAPSE  FOR  SPHERICAL  CHARGES 


The  purpose  of  the  work  discussed  in  this  section  was  to 
verify  that  the  calculational  tools  (i.e.f  computer  codes,  equa¬ 
tions  of  state,  initial  condition,  etc.  )  could  calculate  the  case 
of  the  detonation  of  a  spherical  charge  under  water.  Once  the 
tools  were  validated,  they  could  be  used  to  compute  the  much  more 
difficult  problem  of  the  detonation  of  a  tapered  charge  under 
water.  Section  4.1  reviews  briefly  what  is  known  about  explosion 
bubbles;  Section  4.2  shows  that  the  computed  results  for  the 
spherical  case*  agree  with  empirical  formulas;  Section  4.3 
discusses  steam  generation,  and  Section  4.4  discusses  what  is 
known  about  multiple  bubble  pulsation. 

4.1  LITERATURE  REVIEW 

Experimental  and  theoretical  research  into  the  dynamics  of 
explosion  bubble  growth  and  collapse  began  during  World  War  I, 
reached  a  peak  during  World  War  II,  and  continued  for  10  to  15 
years  after  the  end  of  World  War  II.  It  was  found  that  upon 
detonation  of  a  spherical  or  near-spherical  charge  in  deep  water, 
the  initial  shock  wave  was  followed  by  a  series  of  secondary 
shock  waves.  Theoretical  treatments  showed  that  the  secondary 
shocks  could  be  traced  to  a  pulsating  bubble  that  occupied  the 


‘Detonation  of  a  454-kg  (1000-lb)  spherical  TNT  charge  in  water 
at  a  depth  of  305  m  (1000  ft). 
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region  surrounding  the  explosive  charge.  Laboratory-scale 
experiments  confirmed  this.  An  excellent  summary  of  early  work 
is  given  in  Reference  13,  which  is  a  compendium  of  British  and 
American  reports  concerned  with  explosion  bubble  dynamics. 

Upon  detonation  of  a  spherical  explosive  charge  in  deep 
water,  a  shock  wave  is  transmitted  to  the  water  and  propagates 
away  from  the  source.  After  passage  of  the  initial  shock,  the 
water  retains  a  residual  velocity.  The  detonation  products  ex¬ 
pand  into  the  approximately  spherical  void  around  the  charge;  the 
void  is  caused  by  the  radial  motion  of  the  water  away  from  the 
charge.  Hydrostatic  pressure  halts  the  growth  of  the  explosion 
bubble  eventually,  but  the  pressure  within  the  bubble  at  this 
first  bubble  maximum  radius  (Am)  is  far  less  than  the  hydrostatic 
pressure.  This  pressure  imbalance  is  remedied  through  the  col¬ 
lapse  of  the  bubble,  which  recompresses  the  explosive  products. 
The  pressures  reached  within  the  bubble  at  the  time  of  the  first 
bubble  minimum  radius  (T^)  exceed  the  hydrostatic  pressure;  this 
leads  to  a  second  shock  propagation  radially  outward  in  the  water 
as  the  bubble  reexpands.  This  second  shock  is  the  bubble  pulse. 
The  process  of  bubble  expansion  and  contraction  can  be  repeated 
many  times.  Figure  4.1  (Reference  14)  shows  the  bubble  radius 
versus  time  for  a  specific  case. 

Theoretical  models  of  bubble  oscillation  generally  treat  the 
water  as  an  incompressible  fluid.  Cole  (Reference  14)  and  Snay 
(Reference  15),  among  others,  describe  models  that  ignore  the  gas 
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pressure  within  the  bubble.  Pritchett  (Reference  16)  reviews 
similar  models  that  include  the  bubble  pressure,  and  gives  an 
excellent  summary  of  spherical  bubble  dynamics.  In  addition  to 
the  above  analytic  techniques,  hydrodynamic  calculations  have 
been  performed  using  incompressible  flow  computer  codes  (Refer¬ 
ences  16-18).  The  incompressible  approach  cannot  give  correct 
results  for  the  initial  bubble  expansion,  or  near  the  bubble 
minimum,  because  the  pressures  during  those  times  are  high  enough 
that  water  is  slightly  compressible  (Reference  15). 

Formulas  derived  from  the  analytic  models  can  be  used  to 
predict  Am  and  T^  These  are  found  to  be  a  function  of  charge 
weight,  W,  the  type  of  explosive  used,  and  the  hydrostatic 
pressure,  PQ.  PQ  can  be  expressed  in  terms  of  Z,  the  charge 
depth  plus  the  atmospheric  head: 


A 


m 


W 


1/3 


,1/3 


(4.1) 


and 


T 


1 


W 


1/3 


,5/6 


(4.2) 


The  constants  J  and  K  depend  on  the  units  chosen,  and  on  the 
type  of  explosive  used.  In  m-kg-s  units, 
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(4.3) 


J 


3.50  m4/3/kg 


1/3 


and 


K  =  2.11  s-m5/6/kg1/3  (4.4) 

for  TNT, and  Z  equals  the  water  depth  plus  10  m.  Thus,  for  a 
454-kg  (1000-lb)  spherical  TNT  charge  detonated  at  a  depth  of 
305  m  (1000  ft),  Am  =  3.95  m  (13  ft)  and  Tx  =  134  ms.  For  HBX-1, 

J  =  3.95  m4/3/kg1/3  (4.5) 


and 


v  o  5/6  „  1/3 

K  =  2.41  s-m  /kg  , 


(4.6) 


so  for  a  454-kg  (1000-lb)  spherical  HBX-1  charge  detonated  at  the 
same  depth,  Am  =  4.46  m  (14.6  ft)  and  T^  =  153  ms.  The  above, 
and  many  other  empirical  formulas,  are  summarized  in  a  recent 
report  by  Swisdak  (Reference  19). 

Thus,  the  dynamics  of  explosion  bubble  growth  and  collapse 
in  deep  water  (i.e.,  where  the  hydrostatic  gradient  over  the 
maximum  bubble  diameter  is  small)  are  treatable  in  one  dimension 
using  analytical  tools,  and  these  give  results  that  agree  very 
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well  with  experimental  data  over  the  first  bubble  period  (Refer¬ 
ence  20).  Theoretical  work  concerning  nonspherical  explosion 
bubble  expansion  and  collapse  has  been  limited,  although  the 
analytic  equations  have  been  formulated  (Reference  14).  The  main 
reason  is  that  these  equations  are  extremely  difficult  to  solve 
for  specific  cases. 

4.2  ONE-DIMENSIONAL  VALIDATION  CALCULATIONS 

One-dimensional  (ID)  computer  calculations  were  performed 
simulating  the  detonation  of  454-kg  (1000-lb)  spherical  charges 
of  TNT  and  HBX-1  explosives  at  a  depth  of  305  m  (1000  ft)  in  deep 
water.  Figure  4.2  shows  the  zoning  and  initial  conditions  for 
these  calculations.  The  zoning  is  the  result  of  a  series  of 
preliminary  ID  calculations  in  which  finer  zoning  was  used,  both 
in  the  explosive  and  in  the  water. 

Results  of  these  calculations  include  the  pressure  and  par¬ 
ticle  velocity  at  and  behind  the  initial  shock  in  the  water  at 
all  ranges  as  a  function  of  time,  the  pressure  of  the  explosive 
products  (within  the  explosion  bubble)  versus  time,  and  the 
motion  of  the  water/explosion  products  interface.  For  this 
effort,  which  concentrated  on  the  behavior  of  the  explosion 
bubble,  the  latter  two  results  were  of  the  greatest  interest. 
Figures  4.3a  and  4.3b  show  the  computed  explosion  bubble  radius 
and  bubble  pressure  histories  for  the  spherical  TNT  charge.  The 
bubble  radius  reaches  a  maximum  of  4.05  m  (13.3  ft),  and  then 
decreases  to  a  minimum  radius  of  about  0.9  m  (  3  ft)  at  T^  =145 
ms.  The  internal  bubble  pressure  decreases  rapidly  to  a  value 
much  less  than  the  hydrostatic  pressure  Pc.  The  minimum  bubble 
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3.1  MPa 


pressure  is  about  0.2  MPa  (29  psi),  whereas  PQ  = 

(456  psi).  Bubble  collapse  occurs  because  PQ  reverses  the  radial 
velocity  field  in  the  water.  Since  the  pressure  within  the 
bubble  is  much  less  than  PQ,  it  initially  offers  very  little 
resistance  to  the  collapse  process.  As  the  bubble  is  recom¬ 
pressed  during  collapse,  the  internal  pressure  rises,  and  reaches 
a  value  of  about  30  MPa  (4400  psi)  at  T-^ . 

Table  4.1  compares  computed  and  empirical  (Reference  19) 
values  of  Am  and  T^  for  TNT  and  HBX-1.  For  TNT,  both  A^  and  T^ 
agree  very  well  with  empirical  formulas.  This  is  not  the  case 
for  HBX-1,  however,  where  both  Aj^  and  T^  are  much  larger  than 
empirical  formulas  predict.  After  a  detailed  investigation,  we 
concluded  that  current  knowledge  about  the  HBX-1  EOS  is  not 
adequate  to  enable  calculation  of  explosion  bubble  growth  and 
collapse.  Since  the  computer  simulations  were  validated  for  TNT, 
we  concluded  that  all  2D  calculations  should  be  performed  using 
TNT. 


Similar  calculations  of  explosion  bubble  growth  and  collapse 
have  been  performed  by  Mader  (Reference  21)  and  Walker  (Refer¬ 
ence  22).  For  spherical  Pentolite  charges,  Walker  investigated 
explosion  growth  and  collapse  at  water  depths  of  4267  m 
(14,000  ft)  and  7010  m  (23,000  ft)  and  found  good  agreement  with 
empirical  formulas.  This  work  tends  to  support  the  contention 
that  HBX-1  is  not  characterized  well  enough  for  this  type  of 
computational  effort. 

Computed  results  of  the  shock  in  the  water  were  also  exam¬ 
ined  in  the  ID  calculations.  Figure  4.4  shows  computed  pressure 
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TIME,  ms 


Figure  4.4  Calculated  shock  wave  ra  water  for  a  454-kg 
(1000-lb)  spherical  TNT  charge  detonated  at 
305  m  (1000  ft)  depth  in  water. 


versus  time,  15-m  (49-ft)  and  30-m  (98-ft)  from  a  454-kg 
(1000-lb)  spherical  TNT  charge.  The  initial  shock  is  followed  by 
a  "negative  phase,"  where  the  water  pressure  is  less  than  P0,  and 
this  is  followed  by  the  bubble  pulse.  Figure  4.5  shows  the 
general  characteristics  of  the  pressure  pulse  from  a  deep 
explosion,  and  defines  terms.  These  characteristics  can  be 
predicted  using  empirical  formulas  (Reference  19),  and  the 
computed  results  are  compared  to  the  empirically  predicted  ones 
in  Table  4.2.  Computed  maximum  pressures  are  uniformly  low 
compared  with  empirically  predicted  maximum  pressures.  This  is 
to  be  expected  because  a  finite-difference  computer  code  will 
spread  out  the  shock.  Sternberg  and  Walker  (Reference  23)  found 
that  a  shock-following  technique  was  required  to  accurately 
compute  maximum  shock  pressures  (Ppp)  for  an  underwater 
detonation. 

4.3  THE  PROBLEM  OF  STEAM  GENERATION 

The  water  EOS  of  Walker  and  Sternberg  does  not  compute  the 
production  of  steam  below  the  saturation  curve.  In  their  paper, 
Walker  and  Sternberg  (Reference  11)  examine  the  contribution  of 
steam-generated  pressure  within  the  bubble  for  Pentolite  explo¬ 
sive  detonated  in  water.  Their  conclusion  is  that  steam,  even  if 
generated,  would  not  significantly  raise  the  pressure  within  the 
bubble.  For  our  tapered  charge  calculations,  we  decided  to 
further  investigate  the  problem  of  steam  generation.  A  cylindri- 
cally  symmetric  ID  calculation  was  performed  to  examine  the  peak 
pressures  close  to  the  cylindrical  charge.  In  the  calculation, 
the  explosive  was  detonated  in  over  10  zones,  and  the  first  water 
zone  was  mass-matched  to  the  mass  of  the  last  explosive  zone. 
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variables  are  independent  of  range. 


The  charge  radius  was  set  at  14.75  cm  (corresponding  to  the 
thickest  part  of  the  cylindrical  tapered  charge),  the  first  water 
zone  was  2  cm  thick,  and  an  increasing  zoning  ratio  of  2  percent 
was  used  from  there  on  within  the  water. 

Figure  4.6  plots  the  maximum  shock  pressure  in  the  water 
versus  range  from  the  charge  from  this  calculation.  Only  three 
zones  show  a  maximum  pressure  of  greater  than  5  GPa  (50  kbar); 
for  less  than  that  value,  steam  will  not  form,  even  in  an  equili¬ 
brium  situation,  according  to  the  Butkovich  water  model.  A  cy¬ 
lindrical  envelope  surrounding  the  charge  with  an  approximate 
thickness  of  6  cm  could  possibly  generate  steam  when  the  pressure 
within  that  region  of  shocked  water  falls  below  the  corresponding 
pressure  on  the  saturation  curve.  This  corresponds  to  a  pressure 
of  approximately  1  MPa  (10  bars),  or  about  one-third  of  the  hy¬ 
drostatic  pressure  at  a  depth  of  305  m.  Thus,  steam  production 
can  contribute  (if  it  forms  at  all)  only  to  the  late-time  pres¬ 
sure  within  the  bubble,  i.e.,  near  the  time  when  the  bubble 
reaches  its  maximum  radius.  This  means  that  it  will  not  contri¬ 
bute  significantly  to  the  initial  bubble  growth.  For  these  cal¬ 
culations,  this  is  a  conservative  assumption. 

This  calculation  was  continued  to  investigate  explosion 
bubble  growth  and  collapse  in  cylindrical  symmetry.  This  corres¬ 
ponds  to  an  infinitely  long  cylindrical  charge  with  a  diameter  of 
0.30  m  (1.0  ft).  Figures  4.7a  and  4.7b  show  the  bubble  radius 
and  pressure  versus  time.  If  the  baseline  tapered  charge  were 
long  enough  that  cylindrical  symmetry  dominated,  one  would  expect 
these  results  to  be  valid.  Although  Amax  is  smaller  than  for  the 
spherical  charge  (3.15  m  versus  4.05  m),  T-^  is  much  longer 
(250  ms  versus  145  ms). 
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Figure  4. 


6  Maximum  pressure  versus  range  from  a 
cylindrical  TNT  charqe  with  an  initial 
radius  of  14.75  cm  [  the  maximum  radius 
of  the  chosen  454-kg  (1000-lb)  tapered 
charge ] . 
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Figure  4.7  Calculated  bubble  parameters  versus  time  for  an  infinitely 

long  cylindrical  TNT  charge  with  a  diameter  of  0.30  m  (1.0  ft) 
detonated  in  water  at  a  depch  of  305  m  (1000  ft) . 


4.4  MULTIPLE  BUBBLE  PULSATION 


Explosion  bubble  reexpansion  after  T^  causes  the  bubble 
pulse  to  propagate  in  the  surrounding  water.  This  process  re¬ 
moves  internal  energy  from  the  bubble.  Other  loss  mechanisms, 
such  as  heat  conduction  and  turbulence,  also  tend  to  remove 
energy  from  the  bubble  by  heating  the  surrounding  water. 

Pritchett  (Reference  16)  provides  an  excellent  discussion  of 
these  energy  loss  mechanisms. 

In  underwater  high  explosive  detonations  at  deep  depths, 
bubble  oscillation  can  be  followed  through  up  to  seven  cycles 
before  the  bubble  is  reduced  to  a  benign  gas  pocket,  or  to  many 
small  bubbles.  In  underwater  nuclear  detonations  at  deep  depths, 
bubble  reloading  does  not  appear  to  be  as  effective  because  the 
energy  loss  per  cycle  from  the  bubble  appears  to  be  larger. 

Figure  4.8  qualitatively  compares  high  explosive  and  nuclear 
bubble  pulsation  (Reference  16).  The  key  difference  between  the 
two  is  that  the  nuclear-generated  explosion  bubble  is  filled 
primarily  with  steam  (generated  by  the  initial  shock  wave),  while 
the  high  explosive  bubble  is  filled  primarily  with  noncondensible 
gaseous  explosive  products.  This  point  is  pursued  later  in  this 
report . 

Since  the  purpose  of  this  work  was  to  examine  explosion 
bubble  pulsation  from  high  explosive  charges,  an  effort  was  made 
to  examine  multiple  bubble  pulsation  in  order  to  estimate  how 
much  the  bubble  is  degraded  by  each  oscillation.  Experimental 
data  that  give  bubble  radius  versus  time  for  more  than  one  bubble 
period  appear  to  be  sparse.  The  best  measurements,  and  the  ones 
most  often  quoted,  are  derived  from  Swift  and  Decius  (Refer¬ 
ence  24).  Figure  4.9  shows  bubble  radius  versus  time  for  a 


0  TIME  AFTER  BURST 


Figure  4.8  The  effects  of  bubble  energy  loss — the 
nonmigrating  case  (from  J.  Pritchett, 
Reference  16). 


0.3-kg  (0.66-lb)  TNT  charge  detonated  in  water  at  a  depth  of  91  m 
(300  ft).  The  first  bubble  radius  and  period  are  by  far  the 
largest,  and  a  significant  amount  of  the  bubble  energy  appears  to 
be  absorbed  in  the  first  bubble  collapse  and  reexpansion  process. 

Table  4.3  examines  the  TNT  test  data  from  Reference  22  for 
the  first  three  periods  of  bubble  oscillation.  Average  values  of 
J  and  K  (defined  by  equations  4.1  and  4.2)  decrease  with  each 
oscillation.  Figure  4.10  plots  these  average  values  versus 
period  number,  n.  We  are  not  certain  if  the  power-law  decay  of  J 
and  K  with  period  number  is  valid  for  n  >  3,  so  the  extrapola¬ 
tions  to  n  =  10  are  shown  by  dashed  lines.  Based  on  the  valid 
part  of  the  curves  (n  =  1  to  3),  the  effective  bubble  energy 
(Weff)  relative  to  the  first  period  for  a  nonmigrating  bubble  can 
be  calculated  using  cube  root  scaling;  this  is  done  in 
Table  4.3.  Weff  for  the  second  period  is  only  about  40  percent 
of  that  on  the  first  period,  and  on  the  third  period  is  about  20 
percent.  If  the  extrapolation  of  Figure  4.9  is  valid,  Weff  is 
reduced  to  at  most  10  percent  of  its  first  period  value  by  the 
6th  period.  Also,  Weff  based  on  J  and  K  give  different  values, 
possibly  indicating  that  the  partitioning  of  energy  within  the 
bubble  reqion  is  different  on  each  oscillation.  Bubble  period 
data  from  Slifko  (Reference  25)  for  TNT  charges  were  also 
analyzed,  and  the  data  are  presented  in  Table  4.4.  The  data  for 
the  3.6-  and  86-kg  (8-  and  57-lb)  charges  generally  agree  with 
the  data  presented  in  Table  4.3;  the  0.45-kg  (1-lb)  charges 
appear  to  be  less  efficient  than  the  larger  ones. 

Bubble  energy  losses  arise  from  two  sources:  the  bubble 
pulse  in  the  water  due  to  the  compressibility  of  water,  and 
turbulence  and  related  effects  that  operate  near  the  bubble 
minimum.  To  assess  the  telative  effects  of  these  two  loss 
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Period  number,  n 


Figure 


.10  Jn  and  Kn  versus  n  (data 
from  Swift  and  Decius, 
Reference  24) . 
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mechanisms,  the  ID  spherical  calculation  of  the  detonation  of  a 
454-kg  (1000-lb)  spherical  TNT  charge  in  water  at  a  depth  of 
305  m  (1000  ft)  discussed  in  Section  4.1  was  continued  through 
the  second  bubble  period.  The  results  are  given  in  Figure  4.11, 
which  shows  the  computed  bubble  radius  versus  time.  Calculation 
of  We££  for  T2  relative  to  T^  shows  that  the  relative  bubble 
energy  during  the  second  period  is  60  to  70  percent  of  that  on 
the  first  bubble  period.  This,  then,  is  the  adiabatic  energy 
loss  due  to  the  creation  of  the  bubble  pulse,  since  turbulence 
and  heat  conduction  effects  were  not  included  in  the  calculation. 
Therefore,  at  least  half  of  the  bubble  energy  loss  appears  to  be 
due  to  turbulence  effects  near  the  bubble  minimum,  and  these  act 
most  effectively  during  the  first  bubble  collapse.  These  results 
agree  with  the  theoretical  analysis  given  by  Snay  (Reference  15). 
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Figure  4.11  Computed  explosion  bubble  radius  versus 
time  through  the  second  bubble  period 
for  a  454-kg  (1000-lb)  spherical  TNT 
charge  detonated  in  water  at  a  depth  of 
305  m  (1000  ft) . 


50 


SECTION  5 


TWO-DIMENSIONAL  CALCULATION  OF  THE  EXPLOSION  BUBBLE 
GROWTH  AND  COLLAPSE  CAUSED  BY  THE  DEEP  UNDERWATER 
DETONATION  OF  THE  BASELINE  TAPERED  CHARGE 


A  2D  calculation  was  performed  for  the  baseline  454-kg 
(1000-lb)  TNT  tapered  charge  (Figure  2.4)  using  the  PISCES  2DELK 
(Reference  26)  continuum  mechanics  finite  difference  computer 
code.  The  basic  geometry  for  the  calculation  is  shown  in  Fig¬ 
ure  5.1.  Section  5.1  discusses  the  primary  assumptions  made, 
Section  5.2  presents  the  results,  and  Section  5.3  compares  the 
computed  average  bubble  parameters  with  those  of  an  equal  weight 
spherical  TNT  explosive  charge. 

5.1  CALCULATIONAL  ASSUMPTIONS 


The  following  assumptions  were  made  in  order  to  perform  the 
2D  calculation  in  cylindrical  symmetery  with  the  tapered  charge 
parallel  to  the  sea  surface: 


1.  The  steel  case  enclosing  the  explosive  was  neglected. 
This  may  have  an  effect  on  the  initial  water  shock  pro- 
duced  by  the  charge,  but  should  be  relatively  unimpor¬ 
tant  for  the  explosion  bubble  growth  and  collapse 
process . 

2.  The  hydrostatic  pressure,  P^  was  assumed  constant  over 
the  vertical  maximum  dimensions  of  the  bubble.  The 
worst  case  would  be  if  the  tapered  charge  bubble  had  the 
same  maximum  vertical  dimension  as  an  equal  weight 
spherical  charge.  From  Section  4.1,  Ani  for  a  454-kg 
(1000-lb)  spherical  charge  detonated  at  305  m  (1000  ft) 
is  about  4  m  (13  ft).  PQ  at  305  m  is  about  3.1  MPa 
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Figure  5.1  Geometry  for  2D  calculations. 


(456  psi);  the  variation  over  the  maximum  bubble 
dimensions  is  0.08  MPa  (12  psi),  or  about  2.5  percent  of 
PQ.  Therefore,  the  gradient  in  the  hydrostatic  head  can 
be  neglected.  Inclusion  of  this  gradient  would  make  the 
calculation  fully  three-dimensional. 

3.  Gravity  was  assumed  to  have  little  effect  on  the 

results,  other  than  to  produce  the  hydrostatic  head. 

This  essentially  assumes  that  over  the  first  bubble 
period  the  bubble  is  nonmigrating  (Snay,  Refer¬ 
ence  15).  Using  analytic  formulas,  one  can  estimate  the 
total  amount  of  bubble  migration,  AZ,  for  a  bubble 
caused  by  the  detonation  of  454  kg  (1000-lb)  of  HBX-1 
explosive  at  a  depth  of  305  m  (1000  ft).  This  is  about 
0.9  m  (3  ft),  compared  to  an  Am  of  4.4  m  (14.5  ft). 
Therefore,  the  amount  of  bubble  migration  over  T^  is 
only  21  percent  of  the  maximum  bubble  radius. 

Figure  5.2  shows  the  initial  zoning  and  material  boundaries 
for  the  calculation.  The  coupled  Eulerian-Lagrangian  logic  of 
PISCES  2DELK  was  used — the  Eulerian  grid  described  the  high  ex¬ 
plosive  of  the  tapered  charge  and  the  expansion  cavity  (the 
bubble);  the  Lagrangian  grid  described  the  water.  The  water  had 
a  constant  hydrostatic  pressure  of  3.1  MPa  (31  bars),  consistent 
with  a  water  depth  of  305  m  (1000  ft).  Most  of  the  Eulerian  grid 
was  initially  covered  by  the  Lagrangian  grid  as  Figure  5.2  shows, 
and  the  tapered  charge  position  is  indicated.  The  total  mass  of 
the  tapered  charge  was  449  kg  (990  lb),  and  the  total  length  was 
5.3m  (17. 5  ft) . 

5.2  RESULTS 


The  tapered  charge  was  detonated  from  left  to  right,  as 
indicated  in  Figure  5.2,  and  the  detonation  was  completed  at  a 
time  of  about  0.8  ms.  The  shock  wave  was  rapidly  transmitted  to 
the  water,  and  the  explosion  bubble  began  to  grow.  Bubble  growth 


53 


to  the  time  of  the  first  bubble  maximum  (which  occurs  at  about 
75  ms)  is  shown  in  Figures  5.3/  5.4,  and  5.5.  At  75  ms.  Fig¬ 
ure  5.5  shows  that  the  bubble  is  still  somewhat  elliptical,  so 
the  effect  of  the  large  length-to-diameter  ratio  of  the  tapered 
charge  is  retained  by  the  bubble.  After  75  ms,  the  hydrostatic 
pressure  begins  to  collapse  the  bubble,  as  shown  in  Figure  5.6 
and  5.7.  Collapse  to  first  minimum  (Tj)  occurs  at  approximately 
133  ms. 

Zoning  effects  were  investigated  by  performing  another  cal¬ 
culation  with  twice  the  number  of  zones  close  to  the  charge  to  a 
time  of  25  ms.  Comparison  of  the  computed  results  of  the  two 
calculations  at  that  time  revealed  no  essential  differences, 
indicating  that  the  original  zoning  was  capable  of  resolving  the 
initial  bubble  growth  in  sufficient  detail. 

5.3  COMPARISON  OF  THE  TAPERED  CHARGE  AVERAGED  EXPLOSION  BUBBLE 

PARAMETERS  WITH  THOSE  OF  AN  EQUAL  WEIGHT  SPHERICAL  CHARGE 

Figure  5.7  shows  that  at  T^  the  explosion  bubble  is  "almost 
spherical,"  while  near  the  time  of  maximum  growth  of  the  bubble, 
the  shape  of  the  bubble  was  definitely  elliptical.  It  was  of 
interest  to  know  whether  the  averaged  bubble  parameters  were 
closer  to  those  formed  by  a  spherical  charge  detonated  under¬ 
water,  or  to  an  infinitely  long  cylindrical  charge.  We  had 
previously  calculated  both  symmetries  using  ID  codes  (see 
Section  4). 
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Figure  5.3  PISCES  2DE1 


Figure  5.4  PISCES  2DELK  first  tapered  charge  calculation,  t  =  50  ms 
continued  expansion  of  bubble. 
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Figure  5.7  PISCES  2DELK  first  tapered  charge  calculation,  t  =  133  ms 
the  time  of  the  first  bubble  minimum  (T-i  )  . 


Both  the  total  volume  of  the  explosion  bubble  and  the  aver¬ 
age  pressure  within  the  bubble  were  monitored  in  the  calcula¬ 
tion.  From  the  volume,  an  average  bubble  radius  was  calcualted 
under  the  assumption  that  the  bubble  was  spherical.  Figure  5.8 
compares  this  average  spherized  radius  versus  time  with  the 
bubble  radii  from  an  equal  weight  spherical  charge,  and  the  in¬ 
finitely  long  cylindrical  charge  with  equal  charge  weight  per 
unit  length. 

The  tapered  charge  average  bubble  radius  appears  to  follow 
closely  the  behavior  of  the  bubble  for  a  454-kg  (1000-lb) 
spherical  charge,  but  with  a  slightly  smaller  maximum  radius 
(Am  =  3.8  m  instead  of  4.1  m  for  the  spherical  charge)  and  a 
slightly  shorter  bubble  period  (T^  =  133  ms  instead  of  140  ms), 
consistent  with  the  smaller  maximum  radius. 

Figure  5.9  presents  the  averaged  pressure  within  the  explo¬ 
sion  bubble  versus  time.  The  average  pressure  within  the 
explosion  bubble  for  the  tapered  charge  is  much  closer  to  the 
pressure  within  the  explosion  bubble  of  the  spherical  charge  than 
it  is  for  the  cylindrical  charge.  Also,  it  was  seen  in  the  2D 
calculation  that  there  was  some  variation  of  the  pressure  over 
the  axial  dimensions  of  the  bubble  (this  variation  is  shown  in 
the  figure).  The  variation  is  much  larger  at  the  earlier  times, 
i.e.,  during  the  bubble  expansion,  than  it  is  after  the  bubble 
has  reached  its  maximum  radius  and  begun  to  collapse. 

We  concluded  from  this  comparison  that  (at  least  to  first 
order)  the  bubble  from  the  tapered  charge  could  be  treated  as  a 
spherical  charge,  perhaps  with  the  total  energy  of  the  equivalent 
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BUBBLE  RADIUS 


Figure  5.8  Explosion  bubble  radius  versus  time  for  a  454-kg 
(1000-lb)  spherical  TNT  charge  and  cylindrical 
charge  (diameter  =  0.3  m)  with  volume-averaged 
radius  of  the  bubble  produced  by  the  baseline 
454-kg  tapered  charge  (PQ  =  305  m  hydrostatic 
head) . 
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spherical  charge  slightly  reduced  (by  about  20  percent).  This 
assumption  was  useful  for  the  next  part  of  our  effort,  which 
concentrated  on  injecting  propellant  gasses  into  the  explosion 
bubble  in  an  attempt  to  mitigate  the  bubble  collapse  (see 
Section  6). 

The  fact  that  a  tapered  charge  detonation  deep  under  water 
appeared  to  produce  an  explosion  bubble  with  characteristics  very 
similar  to  those  of  an  equal  weight  spherical  charge  leads  one  to 
suspect  that  such  a  detonation  would  produce  at  least  one  bubble 
pulse.  This  very  important  conclusion  is  corroborated  by  exper¬ 
imental  data  (Reference  27)  from  two  tapered  charges  detonated  in 
deep  water. 
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SECTION  6 


INVESTIGATION  OF  PROPELLANT-ASS iSTED 
PREVENTION  OF  EXPLOSION  BUBBLE  COLLAPSE 


At  the  beginning  of  the  program,  we  thought  that  one  way  to 
eliminate  the  bubble  pulse  caused  by  the  deep  underwater  deto¬ 
nation  of  a  tapered  charge  (if  bubble  pulsation  occurred)  would 
be  to  prevent  the  bubble  from  collapsing.  The  calculation  of  the 
baseline  tapered  charge  discussed  in  the  preceding  section  showed 
that  the  explosion  bubble  caused  by  the  detonation  had  character¬ 
istics  very  similar  to  the  bubble  caused  by  the  detonation  of  a 
spherical  charge  of  the  same  explosive  weight.  Therefore,  one 
would  expect  the  bubble  pulses  to  be  similar  also.  To  prevent 
collapse  of  the  bubble,  we  postulated  that  an  additional  propel¬ 
lant  source  might  be  used  to  supply  gas  to  the  bubble.  Recalling 
the  discussion  of  Section  4.1,  the  pressure  within  the  explosion 
bubble  is  well  below  the  hydrostatic  pressure  at  the  time  of 
maximum  bubble  radius  (Figure  4.3b);  this  pressure  differential 
between  the  bubble  interior  and  the  surrounding  water  is  the 
cause  of  bubble  collapse.  The  purpose  of  the  propellant  source 
would  be  to  raise  the  interior  bubble  pressure,  thereby  halting 
this  process. 

We  realized  early  in  the  proyram  that,  because  of  the  size 
of  the  explosion  bubble,  the  amount  of  propellant  required  would 
be  large.  At  maximum  radius,  the  explosion  bubble  volume  is 
about  270  m  (9500  ft).  Previous  work  with  propellants  (Refer¬ 
ence  28)  was  used  to  estimate  how  much  propellant  would  be 
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was  found  that  the  mass  of  propellant  required  to  raise  the  inte¬ 
rior  bubble  pressure  to  the  hydrostatic  pressure  level  was  pro¬ 
hibitively  large.  A  series  of  ID  calculations  was  then  performed 
to  see  if  smaller  amounts  of  propellant  gas  would  substantially 
mitigate  the  bubble  pulse.  These  are  discussed  in  Section  6.2. 


6.1  ESTIMATE  OF  THE  AMOUNT  OF  PROPELLANT  REQUIRED  TO  AVERT 
BUBBLE  COLLAPSE 


For  a  single  propellant  species  burning  in  a  volume,  Vp  the 
pressure,  P(t)  is  given  by: 


P(t)  =  - ft.  F  -  where  (6.1) 

V  -  “  +  Mt>  -  cj 

M  =  mass  of  the  propellant 
F  =  propellant  impetus  (energy  content) 

C  =  covolume  of  the  propellant  gasses 
p  =  initial  density  of  the  propellant 
A( t )  =  fraction  of  the  propellant  burned. 

Setting  A  =  1  at  t  =  «  gives 


P 

00 


MF 

V  -  MC 


(6.2) 


Properties  for  a  typical  U.S.  Navy  gun  propellant,  SPDN  (Refer¬ 
ence  24)  were  used:  F  =  9.35  x  103  J/kg,  and  C  = 

9.33  x  10“4  m3/kg . 
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The  amount  of  propellant  required  to  raise  the  interior 
bubble  pressure  to  the  hydrostatic  pressure,  PQ,  at  the  time  of 
the  first  bubble  maximum  for  the  454-kg  (1000-lb)  charge  can  be 
calculated  as  follows.  First, 


Po  "  Pbubble  +  P» 


(6.3) 


where  P^^ie  =  0.2  MPa  (2  bars)  at  the  time  of  first  bubble 
maximum  (see  Figure  4,3b).  This  defines  P#  as 


P  =  3.1  MPa  -  0.2  MPa  =  2.9  MPa  . 
00 


(6.4) 


Next,  the  mass  of  propellant  can  be  calculated  using  Equa 
tion  6.2: 


M 


VP 

00 

F  +  CP 


00 


_ (270  m3) (2.9  x  IQ6  Pa) _ 

9.35  x  105  J/kg  +  (9.33  x  10-4  m3/kg)(2.9  x  106  Pa) 


=  835  kg  (1840  lb). 


(6.5) 
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Thus,  if  propellant  gas  could  be  injected  instantaneously 
into  the  explosion  bubble  at  this  time,  835  kg  (1840  lb)  would  be 
required  to  raise  the  bubble  pressure  to  the  hydrostatic  pressure 
at  a  depth  of  305  m  (1000  ft).  This  propellant  weight  is  almost 
twice  the  weight  of  the  454-kg  (1000-lb)  explosive  charge. 

This  is  clearly  a  prohibitively  large  amount  of  propellant, 
and  thus  it  did  not  appear  feasible  to  completely  arrest  the 
bubble  collapse  process.  Actually,  this  estimate  is  probably  a 
lower  bound,  since  some  bubble  oscillation  might  still  occur. 

6.2  CALCULATIONS  AND  ANALYSIS 

Since  it  did  not  appear  possible  to  avert  bubble  collapse 
with  any  reasonable  amount  of  propellant,  an  effort  was  made  to 
examine  to  what  extent  a  reasonable  amount  of  propellant  could 
reduce  the  bubble  pulse.  This  required  running  ID  spherical 
calculations  well  past  the  time  of  the  first  bubble  minimum  so 
that  the  bubble  pulse  could  be  examined  at  a  reasonable  range  in 
water.  The  ID  calculations  were  adequate  because  of  the  simi¬ 
larity  of  the  tapered  charge  explosion  bubble  to  a  spherical 
charge,  discussed  in  Section  5.3. 

Calculations  were  performed  wherein  propellant  gasses  were 
injected  into  the  explosion  bubble  at  constant  rates.  Figure  6.1 
shows  the  computational  method  that  we  followed.  In  all  cal¬ 
culations,  a  454-kg  (1000-lb)  spherical  TNT  charge  -as  initially 
detonated.  The  propellant  source  was  placed  within  the  explosion 
bubble  using  a  special  boundary  condition  to  inject  the  propel¬ 
lant  gas  into  the  growing  explosion  bubble.  Table  6.1  summarizes 
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Figure  6.1  Sketch  showing  ID  computational  method 
for  adding  a  constant  propellant  gas 
source  to  the  TNT  explosion  bubble. 
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Table  6.1  Summary  of  initial  conditions  for  the  ID 
calculations  with  propellant  qas  mass 
addition. 


Initial  Conditions  Used  in  All  of  the  Calculations 


TNT:  453-kg  (1000-lb)  charge 

Water:  Constant  hydrostatic  pressure 

Propellant  Source 

Density  of  gasses: 


Exit  velocity 
of  gasses: 

Internal  energy 
of  gasses: 

Pressure  of  gasses: 


3.1  MPa  (460  psi) 

0.1  g/cc 

1900  m/s  (6200  ft/s) 
2570  J/g 

7.7  MPa  (1130  psi) 


Initial 

Conditions  Unique  to  Each 

Calculation 

Calculation 

rQ,  cm 

Mass  Addition  Rate 

kg/ms 

lb/ms 

1 

0.0 

0.0 

(0.0) 

2 

0.573 

0.0784 

(0.173) 

3 

1.14  6 

0.  3136 

(0.691) 

4 

4.584 

5.018 

(11.06) 

5 

9.168 

20.070 

(44. 25) 
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the  initial  conditions  for  these  ID  spherical  calculations.  The 
propellant  gasses  were  given  a  density,  an  exit  velocity,  an  in¬ 
ternal  energy,  and  a  pressure  typical  of  those  quantities  for 
operating  rockets  (Reference  29).  The  mass  injection  rate  was 
varied  by  varying  the  radius  rQ,  as  indicated  in  Figure  6.1,  over 
which  the  propellant  source  was  injected  into  the  explosion 
bubble.  The  mass  injection  rates  ranged  from  0.0  kg/ms  (no  mass 
addition  at  all)  to  20.07  kg/ms,  as  shown  in  Table  6.1.  We  felt 
that  this  series  of  calculations  would  give  us  some  indication  of 
whether  the  addition  of  hot  propellant  gasses  to  the  explosion 
bubble  would  contribute  to  mitigation  of  the  bubble  pulse. 

Results  of  this  series  of  calculations  are  given  in  Fig¬ 
ure  6.2,  which  shows  the  computed  bubble  radius  versus  time. 

Some  "cushioning"  of  the  explosion  bubble  is  seen,  even  with  the 
lowest  (0.0784  kg/ms)  propellant  injection  rate.  Too  large  a 
rate  (20.07  kg/ms)  leads  to  the  formation  of  a  much  larger 
bubble,  and  this  probably  not  useful.  An  injection  rate  of  about 
5  kg/ms  appears  to  minimize  bubble  collapse,  while  not  signifi¬ 
cantly  enlarging  the  bubble  maximum  radius.  This  "optimized" 
injection  rate  leads  to  the  addition  of  an  even  larger  mass  of 
propellant  to  the  explosion  bubble,  about  1,000  kg  (2200  lb)  over 
one  bubble  period  (220  ms).  The  constant  mass  injection  rate  is 
the  most  feasible  method  of  supplying  hot  gasses  to  the  explosion 
bubble;  clearly,  the  "instantaneous"  method  of  Section  6.1  is 
unrealistic.  The  gas  is  then  allowed  to  contribute  to  bubble 
growth,  however,  thus  lengthening  the  bubble  period  and  requiring 
even  more  propellant. 
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BUBBLE  RADIUS, 


(1000-lb)  spherical  TNT  charge  detonated  at  a 
depth  of  305  m  (1000  ft)  in  water;  compared  for 
various  constant  injection  rates  of  propellant 
gasses  into  the  explosion  bubble. 


Figure  6.3  shows  pressure  versus  time  at  a  range  of  40  m  for 
the  case  of  no  propellant  gas  injection,  and  the  "optimized"  case 
where  5.0  kg/ms  (11.1  lb/ms)  of  propellant  gas  is  injected  into 
the  explosion  bubble.  The  initial  shock  wave  is  not  affected  by 
the  propellant  injection,  but  the  bubble  pulse  peak  overpressure 
is  reduced  from  0.7  MPa  (105  psi)  to  0.3  MPa  (45  psi),  more  than 
a  factor  of  two,  by  the  propellant  gas  injection.  Total  positive 
phase  impulse  in  the  bubble  pulse  is  decreased  22  percent  from 
22.2  MPa-ms  (3330  psi-ms)  to  17.4  MPa-ms  (2600  psi-ms). 

These  results  were  encouraging,  but  they  also  showed  that 
even  to  mitigate  the  bubble  pulse  is  very  difficult.  Injection 
of  propellants  into  the  explosion  bubble  appears  to  cushion 
bubble  collapse,  thereby  reducing  the  peak  bubble  pulse  pressure 
in  the  surrounding  water,  but  at  the  same  time  the  bubble  pulse 
positive  phase  is  spread  out  over  a  longer  period  of  time. 

We  concluded  that  propellant-assisted  bubble  pulse  miti¬ 
gation  was  not  effective  enough  in  reducing  impulse,  considering 
the  large  amount  of  propellant  that  would  somehow  have  to  be 
injected  into  the  bubble  during  its  expansion  and  collapse. 
Further,  the  discussion  of  the  phenomenology  of  multiple  bubble 
pulsation  (Section  4.4)  implies  that  by  cushioning  the  first 
bubble  collapse,  one  might  actually  increase  the  overall 
stability  of  the  bubble.  Thus,  while  decreasing  the  bubble  pulse 
impulse  per  bubble  period,  one  may  actually  increase  the  total 
number  of  bubble  oscillations. 
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Figure  6.3  Calculated  pressure  versus  time  at  the  40-m  (131-ft) 
range  from  a  453-kg  (1000-lb)  spherical  TNT  charge 
detonated  in  water  at  a  depth  of  305  m  (1000  ft) 
compared  to  that  calculated  for  the  TNT  charge  with 
a  constant  propellant  gas  injection  source  inside 
the  explosion  bubble. 
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SECTION  7 


INVESTIGATION  OF  TWO  POSSIBLE  METHODS  OF  MITIGATION 
OF  TAPERED  CHARGE  EXPLOSION  BUBBLE  PULSATION 


We  investigated  two  additional  methods  to  see  if  bubble 
pulsation  could  be  significantly  reduced,  even  though  it  cannot 
be  eliminated  entirely.  Both  methods  essentially  employed  the 
same  principle  of  using  the  first  bubble  collapse  to  destroy  the 
continuity  of  the  explosion  bubble.  The  first  method  was  to  make 
the  explosion  bubble  from  a  conventional  tapered  charge  suffi¬ 
ciently  asymmetric  at  that  turbulence  would  absorb  more 
energy,  leaving  less  for  later  pulsations.  This  is  discussed  in 
Section  7.1,  along  with  an  initial  2D  calculation  for  a  specific 
geometry. 

The  second  method  was  to  fill  the  explosion  bubble  with 
explosion  products  that  are  either  condensible  gasses  (e.g., 
steam)  or  are  solids  soluble  in  water,  by  using  a  different  ex¬ 
plosive  in  the  tapered  charge.  At  T^f  these  products  would 
condense  and/or  dissolve  into  the  water,  preventing  the  rapid 
pressure  rise  that  occurs  with  conventional  charges,  and  thereby 
reducing  bubble  pulsation.  Such  explosives  were  developed  about 
10  years  ago  at  the  U.S.  Naval  Ordnance  Laboratory,  White  Oak, 

MD.  Section  7.2  discusses  these  explosives,  and  how  they  could 
be  used  to  reduce  bubble  pulsation  for  simulated  underwater 
nuclear  detonations. 
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7.1  ASYMMETRIC  BUBBLE  COLLAPSE 


Examination  of  explosively  generated  explosion  bubbles  (Sec¬ 
tion  4)  revealed  that  the  first  bubble  collapse  and  reexpansion 
removed  the  greatest  amount  of  energy  from  the  bubble  (approxi¬ 
mately  60  percent),  and  that  turbulence  and  heat  conduction 
effects  were  responsible  for  more  than  half  of  this  energy 
loss.  This  is  true  for  detonations  in  deep  water:  a  non¬ 
migrating  bubble  is  remarkably  spherical  over  many  periods  of 
oscillation;  i.e.,  it  is  dynamically  stable.  If  the  symmetry  of 
the  explosion  bubble  could  be  destroyed,  particularly  at  T^ , 
instabilities  might  be  enhanced,  thereby  making  the  energy  loss 
mechanisms  even  more  effective.  This  in  turn  would  damp  out 
bubble  pulsation  more  rapidly. 

An  asymmetric  bubble  can  be  initiated  by  placing  an  air- 
filled  "get  lost"  pipe  directly  behind  the  tapered  charge,  as 
shown  in  Figure  7.1.  The  initial  water  shock  at  the  model  loca¬ 
tion  would  not  be  affected,  but  explosive  gasses  will  propagate 
down  the  air-filled  pipe  once  the  charge  has  completely  deto¬ 
nated.  Bubble  growth  would  then  proceed  in  an  asymmetric  manner, 
as  illustrated  in  Figure  7.2.  For  this  technique  to  work,  the 
bubble  must  remain  asymmetric  through  the  collapse  phase  to  the 
time  of  first  minimum  radius  (Tj).  This  is  illustrated  in  Fig¬ 
ures  7.3a  and  b. 

A  2D  calculation  was  performed  to  further  investigate  this 
method  of  bubble  pulse  mitigation.  We  recognized  that  the  entire 
concept  could  not  be  investigated  because  heat  conduction  and 
turbulence  models  were  not  included  in  the  computational 
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Figure  7.1  Tapered  charge  with  pipe  attached  to  introduce 
bubble  asymmetry. 
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P0  =  305  m  hydrostatic  head 
or  “deep  burst," non¬ 
migrating  bubble. 


Figure  7.2  What  asymmetric  bubble  might  look  like  at  first 
maximum. 
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SYMMETRY 

AXIS 

Bubble  collapse  begins  to  occur,  and  Large  velocity  imbalance  occurs, 

bubble  gasses  escape  to  pocket.  creating  extreme  turbulence  near 

bubble  minimum.  This  turbulence 
destroys  the  bubble  by  breaking 
it  into  smaller  bubbles.  Cooling 
occurs  quickly,  and  little  bubble 
reexpansion  occurs. 


Figure  7.3  How  bubble  collapse  might  take  place  using  this 
technique. 


treatment.  Therefore,  the  bubble  collapse  to  could  only  be 
evaluated  qualitatively  to  determine  what  effect  the  bubble 
asymmetry  might  have  on  the  later  oscillations  of  the  bubble. 

Figure  7.4  shows  the  initial  zoning  for  this  second  PISCES 
2DELK  calculation.  The  air-filled  pipe  was  assigned  an  initial 
diameter  of  0.3  m  (1  ft),  the  same  as  the  maximum  diameter  of  the 
tapered  charge.  The  pipe  itself  was  not  modeled  in  the  calcula¬ 
tion,  and  the  open  region  on  the  inside  was  treated  as  a  void. 
Figure  7.5  shows  the  explosion  bubble  and  the  velocity  field  in 
the  surrounding  water  at  5  ms.  Explosive  gasses  had  reached  the 
end  of  the  tube  and  stagnated,  forming  a  second,  smaller  bubble 
and  a  "stagnation  bubble."  By  10  ms  (Figure  7.6),  the  pipe 
region  between  the  explosion  bubble  and  the  stagnation  bubble  hau 
closed  off  due  to  the  higher  pressure  in  the  surrounding  water. 

At  25  ms  (Figure  7.7),  the  stagnation  bubble  had  ceased  to 
expand.  At  that  time,  it  contained  about  one  percent  of  the  TNT 
explosive  products;  the  rest  were  sealed  in  the  tapered  charge 
explosion  bubble. 

Comparison  of  Figure  7.7  with  Figure  5.3  shows  that  the 
explosion  bubble  was  not  significantly  affected  by  the  presence 
of  the  open  pipe.  The  velocity  vectors  of  Figure  7.7  show  that 
the  velocity  field  in  the  water  is  still  roughly  symmetric  on 
each  end  of  the  tapered  charge  explosion  bubble. 

The  above  results  led  us  to  the  conclusion  that  the  pipe  was 
too  small  in  this  calculation  to  significantly  affect  the  growth 
of  the  explosion  bubble.  Therefore,  it  was  terminated  at  25  ms, 
and  a  third  calculation  was  generated.  This  calculation  was 


80 


tt’21 108 


Figure  7-4  Initial  zoning  of  the  second  tapered  charge  calculation  showing  the 
locations  of  the  charge  and  the  air-filled  pipe. 


calculation  at  10 


similar,  but  the  initial  diameter  of  the  air-filled  pipe  was 
doubled,  to  0.6  m  (2  ft),  and  then  gradually  tapered  to  a  0.3-m 
( 1— f t )  diameter.  The  length  of  the  pipe  was  the  same  as  in  the 
second  calculation.  Figure  7.8  shows  the  geometry  and  zoning  for 
this  calculation.  The  intended  purpose  of  the  larger-diameter, 
air-filled  pipe  was  to  aid  in  generating  a  larger  stagnation 
bubble,  closer  to  the  tapered  charge  explosion  bubble. 

Figures  7.9,  7.10,  and  7.11  are  computer-generated  zoning, 
material  boundary,  and  velocity  vector  plots  at  25,  50,  and 
75  ms,  respectively.  A  much  larger  stagnation  bubble  was  formed 
(Figure  7.9)  and  was  continuing  to  grow  at  25  ms.  It  remained 
separated  from  the  explosion  bubble  by  a  thin  region  of  water, 
and  contained  about  six  percent  of  the  explosive  mass,  clearly  an 
improvement  over  the  results  of  the  second  calculation.  Fig¬ 
ure  7.11  shows  that  the  velocity  field  in  the  water  is  definitely 
asymmetric  at  approximately  the  time  of  the  maximum  size  of  the 
explosion  bubble. 

It  was  uncertain  at  this  point  whethtr  or  not  the  asymmetry 
created  was  sufficient,  and  so  the  calculation  was  continued 
through  Tj^  to  the  early  portion  of  the  bubble  reexpansion.  At 
100  ms  (Figure  7.12)  a  vortex  was  formed  in  the  water  between  the 
expansion  and  stagnation  bubbles.  Whether  or  not  the  explosion 
and  stagnation  bubbles  would  truly  be  disconnected  during  the 
collapse  phase  is  somewhat  in  doubt.  For  example,  the  pipe,  or 
at  least  pieces  of  it  (not  modeled  in  this  calculation)  would 
still  be  contained  in  this  region,  which  is  less  than  0.5  m 
thick.  It  would  represent  a  major  discontinuity  there,  causing 
the  growth  of  instabilities. 
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The  results  of  the  inviscid  continuum  calculation  indicate 
that  the  collapse  and  reexpansion  of  the  explosion  bubble  are 
controlled  largely  by  the  explosion  bubble.  Figure  7.13  shows 
the  results  at  125  ms,  and  Figure  7.14  shows  the  results  at 
140  ms,  which  is  essentially  T^.  The  bubble  collapse  point 
appears  to  be  located  off-axis,  indicating  that  the  bubble  might 
assume  a  toroidal  shape  for  a  short  period  of  time.  Bubble  re¬ 
expansion  is  asymmetric,  as  shown  in  Figure  7.15  (158  ms). 

It  was  obvious  that  one  or  two  calculations  were  insuffi¬ 
cient  to  entirely  solve  the  problem  of  tapered  charge  bubble 
pulse  mitigation;  an  additional  effort  will  be  required  to  ad¬ 
dress  the  more  difficult  problems  of  turbulent  mixing.  However, 
the  results  are  encouraging.  It  appears  that  it  is  possible  to 
enhance  instabilities  in  this  manner.  Another  way  to  further 
investigate  this  concept  of  bubble  pulse  mitigation  would  be  to 
perform  a  test  series  using  tapered  charges  in  deep  water,  some 
with  open  pipes  attached,  and  some  without. 

7.2  STEAM- PRODUCING  EXPLOSIVES 

Two  steam-producing  explosives  were  tested  by  the  U.S.  Naval 
Ordnance  Laboratory,  White  Oak,  MD  (References  30,31).  These 
explosives  were  aluminized  lithium  perchlorate  (LiCl04 Ol^O/Al ) 
and  aluminized  hydrogen  peroxide  (H2O2/AI).  The  basic  detonation 
reactions  are  as  follow: 


3Li CIO  .  •  3 II  „0  +  8A1  3LiCl  +  4  Alo0^  +  9II„0 
4  2  2  3  2 


(7.1) 
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Figure  7.14  Material  boundary  and  vector  velocity  plot  from  the  tjhird 
tapered  charge  calculation  at  140  ms,  the  time  of  the 
first  bubble  minimum. 


and 


3H202  +  2A1  -*•  Al203  +  3H20  .  (7.2) 

The  reaction  products  other  than  steam  are  either  solid  or  are 
soluble  in  water;  i.e.,  no  noncondensible  gasses  are  formed.  The 
aluminized  lithium  perchlorate  (Lithanol)  was  found  to  be  by  far 
the  safer  of  the  two  explosives. 

The  purpose  of  the  program  was  to  develop  an  explosive  that 
would  generate  a  steam-filled  bubble,  as  in  a  nuclear  detona¬ 
tion.  Bubble  pulsation  is  substantially  reduced  in  the  nuclear 
case,  as  we  discussed  in  Section  4.4.  Small  charge  tests  were 
performed  to  obtain  values  of  underwater  shock  wave  parameters  in 
comparison  with  Pentolite  (Reference  30).  Larger  charge  tests, 
22.6  kg  and  136  kg  (50  lb  and  300  lb),  were  performed  to  examine 
bubble  pulsation  and  migration  characteristics  (Reference  31). 
This  section  reviews  the  bubble  pulsation  characteristics  of 
Lithanol . 

Four  experiments  were  conducted  in  which  approximately  23-kg 
(50-lb)  spherical  charges  of  Lithanol  were  detonated  in  shallow 
(46-m-deep)  water  at  a  depth  of  about  30.5  m  (100  ft).  Four 
similar  Pentolite  experiments  were  also  conducted.  These  test 
conditions  are  summarized  in  Table  7.1.  Test  results  showed  that 
the  explosion  bubbles  formed  in  these  tests  were  not  strongly 
migrating . 
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able  7.1  Experimental  conditions  for  23-kg  (50-lb) 
spherical  Lithanol  and  Pentolite  charges 
detonated  in  water  at  a  depth  of  30.5  m 
(100  ft) (data  from  Reference  31). 
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Bubble  pulses  were  recorded  for  at  least  four  bubble  periods 
for  both  Lithanol  and  Pentolite;  Figure  7.16  compares  the  first 
three  bubble  pulses,  as  presented  in  Reference  31.  As  Phillips 
and  Willey  point  out,  "The  shape  of  the  first  bubble  pulse  is 
essentially  the  same  for  both  explosives.  The  second  and  third 
pulses  are  considerably  different,  the  Lithanol  showing  many 
large  spikes  as  opposed  to  the  generally  rounded  appearance  of 
the  Pentolite  pulses.  In  actuality,  even  more  spikes  were 
visible  on  the  original  playouts  than  could  be  shown  on  these 
tracings."  The  pressure  spikes  in  the  Lithanol  waveforms  were 
attributed  to  water  jets  impacting  water  on  the  opposite  side  of 
the  bubble,  much  like  the  "water  hammer"  effect  in  the  case  of  a 
migrating  bubble. 

Table  7.2  compares  the  first  four  measured  bubble  periods; 
the  effective  energies  relative  to  the  first  period  bubble  energy 
(the  cubes  of  the  ratios  of  Tn  to  Tj ,  as  in  Section  4.4)  are  also 
compared.  The  Pentolite  experiments  indicate  a  second  period 
bubble  efficiency  of  about  60  percent;  generally,  bottom  and 
surface  effects  appear  to  help  the  bubble  retain  more  energy  than 
it  would  in  deep  water.  Even  so,  the  relative  efficiency  of  the 
Lithanol-produced  bubbles  (23-25  percent)  is  much  less  than  the 
comparative  Pentolite  data,  on  the  TNT  data  in  deep  water  (about 
40  percent,  Table  4.4).  The  effects  of  steam  condensation  are 
even  more  effective  on  the  later  bubble  periods. 

Thus  it  appears  that  if  Lithanol  charges  are  fired  in  deep 
water,  bubble  oscillation  will  be  damped  out  even  more  rapidly 
than  indicated  by  the  above  data.  Reducing  the  effective  bubble 
energy  by  condensing  the  gasses  and  dissolving  the  rest  of  the 
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THIRD  PULSE 

PENTOLITE  (SHOT  PW45)  LITHANOL  (SHOT  PW-48) 


Figure  7.16  Comparison  of  bubble  pulse  resulting  from  a  23-kg 
(50-lb)  spherical  Pentolite  and  Lithanol  charge 
detonation  in  water  at  a  depth  of  30  m  (100  ft) . 
The  gage  depth  is  13  m  (59  ft) . 
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bubble  contents  in  water  will  lead  to  substantially  mitigated 
bubble  pulses.  Lithanol  is  a  safe,  granular  explosive,  and  it 
appears  that  tapered  charges  could  be  designed  to  use  it.  As 
with  asymmetric  bubble  collapse,  further  calculations  and  experi¬ 
ments  will  be  required  to  quantify  this  method  of  bubble  pulse 
mitigation. 


SECTION  8 


SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


Physics  International  Company  participated  in  a  Defense 
Nuclear  Agency-sponsored  effort  to  examine  bubble  pulsation  re¬ 
sulting  from  the  deep  underwater  detonation  of  tapered  high 
explosive  charges.  The  objectives  of  the  effort  were  to: 

1.  Determine  if  the  explosion  bubble  formed  by  the  deep 
underwater  detonation  of  a  tapered  charge  would  produce 
bubble  pulses. 

2.  Investigate  methods  of  mitigating  or  eliminating  these 
bubble  pulses,  should  they  occur. 

During  the  effort,  nonmigrating  bubble  pulsation  from 
spherical  charge  detonations  was  examined  in  detail,  as  much 
information  is  available  from  past  efforts.  This  information 
includes  measurements  of  bubble  pulsation  that  have  resulted  in 
reliable  empirical  formulas  describing  this  behavior  (at  least 
over  the  first  period  of  oscillation),  and  theoretical  efforts 
that  can  predict  the  general  features  of  bubble  expansion  and 
collapse.  The  current  effort  extended  this  work  by  calculating 
the  expansion  and  collapse  of  a  spherical  explosion  bubble  from 
the  time  of  detonation.  The  computed  first  bubble  pulse  compared 
very  well  with  empirical  formulas,  even  though  turbulence,  water 
jetting,  heat  conduction,  and  other  energy  loss  mechanisms  which 
operate  during  the  bubble  collapse  phase  were  not  modeled. 
Successful  computation  of  spherical  explosion  bubbles  made 


the  calculational  method  credible  for  the  more  difficult  tapered 
charge  explosion  bubble  calculations,  and  also  provided  a 
detailed  look  at  bubble  dynamics. 


Major  results  of  the  effort,  as  related  to  the  objectives, 
are  as  follow: 


1.  At  least  one  bubble  pulse  will  occur  if  a  454-kg 
(1000-lb)  tapered  charge  is  detonated  at  a  depth  of 
305  m  (1000  ft).  The  strength  of  this  pulse  should  be 
almost  as  great  as  if  the  charge  had  been  spherical. 

2.  Elimination  or  mitigation  of  the  bubble  pulse  by  injec¬ 
tion  of  propellant  gasses  into  the  bubble  does  not 
appear  feasible  because  the  weight  of  propellant 
required  would  exceed  the  explosive  charge  weight  by  at 
least  a  factor  of  two. 

3.  Energy  loss  mechanisms  operating  during  the  period  of 
the  first  bubble  minimum  absorb  the  greatest  fraction  of 
the  bubble  energy,  even  for  conventional  (i.e.,  non-nuclear) 
charges.  Therefore,  enhancement  of  bubble  asymmetries  and 
turbulence  might  mitigate  further  bubble  pulsation.  An  air- 
filled  pipe  attached  to  the  tapered  charge  is  one  way  to 
encourage  the  creation  of  an  asymmetric  bubble. 

4.  Steam-producing  explosives,  developed  by  the  U.S.  Navy 
to  simulate  the  characteristics  of  the  nuclear  explosion 
bubble,  have  already  been  shown  to  be  effective  in  miti¬ 
gating  bubble  pulses  in  shallow  water.  Analysis  sug¬ 
gests  that  in  deep  water,  use  of  such  explosives  would 
lead  to  even  more  effective  bubble  pulse  mitigation. 


The  following  recommer dations  are  made  as  a  result  of  this 
initial  effort: 


1.  Turbulence  and  heat  conduction  effects  must  be  incor¬ 
porated  into  the  calculational  technique  in  order  to 
accurately  compute  multiple  bubble  pulsation. 
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2.  Deep  water  testing  of  small  scale  (~  100-lb)  tapered 
charges  should  be  performed.  Some  of  these  tests  might 
incorporate  air-filled  pipes;  bubble  pulse  measurement 
will  show  how  well  this  method  of  bubble  pulse  mitiga¬ 
tion  works. 

3.  Lithanol,  a  steam-producing  explosive,  should  be  further 
investigated  for  use  in  tapered  charges. 
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